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The spin orbit interactions (SOI) of light mediated by single scattering from plasmon resonant metal nanopar- 
ticles (nanorods and nanospheres) are investigated using explicit theory based on Jones and Stokes-Muellcr 
polarimetry formalism. The individual SOI effects are analyzed and interpreted via the Mueller matrix-derived, 
polarimetry characteristics, namely, diattenuation d and retardance 8. The results demonstrate that each of the 
contributing SOI effects can be controllably enhanced by exploiting the interference of two neighboring modes 
in plasmonic nanostructures (orthogonal electric dipolar modes in rods or electric dipolar and quadrupolar 
modes in spheres). 
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Spin orbit interaction (SOI) of light has evoked inten- 
sive theoretical and experimental investigations in the 
past few years owing to fundamental interests and po- 
tential applications in nano-optics [1-8]. The SOI phe- 
nomenon signifies interconversion between the spin angu- 
lar momentum (SAM, circular polarization represented 
by polarization helicity a = ±1) and orbital angular 
momentum (OAM, helical phase of light represented 
by topological vortex charge I = 0, ±1, ±2, . . .), which 
can be mediated via a number of processes involving 
light-matter interactions [2-8]. Among these, while SOI 
arising due to propagation through anisotropic media 
is well known and studied over decades [9], that oc- 
curring in isotropic media has only been noticed re- 
cently [2-8]. For example, such intrinsic coupling be- 
tween SAM and OAM has been observed in tight fo- 
cusing of fundamental or higher order Gaussian beams 
[8], high numerical aperture (NA) imaging in isotropic 
homogenous/inhomogeneous media [6] and in scatte- 
ring from micro and nano systems [4,5]. The spin or- 
bit coupling in such cases have accordingly been mod- 
eled using different methods such as Debye Wolf theory 
for focusing [6,10,11], Mie theory for scattering [4,5]. 
A number of interesting and intricate polarization ef- 
fects associated with topological evolution of phase have 
also been observed and explained by SOI and subse- 
quent conservation of total angular momentum (TAM) 
of light [2,4,6,12-14]. Note each of the individual ef- 
fects associated with SOI are manifested as a measur- 
able change in the spatial polarization characteristics 
of light. For conceptual and practical reasons, model- 
ing SOI via the conventional polarimetry characteris- 
tics may thus prove to be useful [15]. In this letter, 
we present a generalized theory based on conventional 
Jones and Stokes-Mueller formalism for analysis, inter- 
pretation and quantification of the SOI effects (medi- 



ated by scattering) via the individual medium polarime- 
try characteristics, namely, diattenuation d (differential 
attenuation of orthogonal polarization states) and re- 
tardance 6 (phase shift between orthogonal polarization 
states) [15,16]. Importantly, we show using explicit the- 
ory and illustrative examples that spin orbit interactions 
by scattering can be significantly enhanced in plasmon 
resonant metal nanoparticles/nanostructures (nanorods 
and nanospheres). Moreover, in such nano-systems, each 
of the contributing SOI effects can be desirably tuned 
(optimized/enhanced) via the diattenuation and the re- 
tardance parameters by changing the wavelength of light 
and controlling the size and shape of the nanoparticles. 

In order to study SOI of light mediated by the scat- 
tering process, let us choose the right handed Cartesian 
coordinate system with the incident light (plane wave) 
propagating in the Z direction and assume that X and 
Y are the two orthogonal axes, representing the polar- 
ization axes in the laboratory reference frame. The scat- 
tered electric field (E s ) components can be related to the 
incident field (E l ) components in the laboratory frame 
by the action of the following transfer function (J) 

E s « T z {-^)T y {-0)S{9)T z (4>)& = JE\ where, 

(E a + E/j cos 20 Ep sin 20 E 1 cos \ 
.Eg sin 20 E a — Ef3 cos 20 i? 7 sin0 1(1) 
-E 1 cos —E y sin E a + EpJ 

where, 9 is the scattering angle and is the azimuthal 
angle. The transformation matrix T z (<p) transforms the 
incident laboratory frame field vector to the scattering 
plane. The inverse transformation matrices T 2 (— 0) and 
T y {—6) transforms the scattered field from the scatte- 
ring co-ordinate (r, 9, 0) to the laboratory co-ordinate 
(X,Y,Z). The matrix S(9) (defined in the scattering 
plane) includes the effect of the scattering process in its 
elements; S^^) - scattered field polarized parallel to the 
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scattering plane(#), and Si (9) - polarized perpendicular 
to the scattering plane (0) [17]. The scattered field de- 
scriptors, E a (9), Ep(6) and E 1 (9) of Eq. 1, are related 
to the elements of S(9) as 

E a = S 2 cos6 + Si;Ep = S 2 cos6> - Si;E^ = -S 2 sm8 

(2) 

Note that the electric field transformation matrix of Eq. 
1 is similar in nature to that for focusing [6,8, 10]. Thus 
akin to focusing, the SOI effect (SAM to OAM con- 
version) can readily be verified by applying the Jones 
vector [1 i 0] T of right circularly polarized (RCP) light 
on the matrix of Eq. 1 and by decomposing the result- 
ing field into three uniform polarization components (as 
previously done for focusing of fundamental Gaussian 
laser beam [10]). However, unlike focusing, the matrix 
S(0) (in Eq. 1) additionally incorporates information on 
the interactions of light with micro or nano scale ob- 
jects and thus the SOI mediated by scattering from such 
systems is expected to be more complex. Since in an 
actual experiment, the set of analyzing optics and de- 
tectors (kept in the X — Y plane, in far field) typically 
detect the transverse components of the scattered fields, 
henceforth we consider the transverse components of the 
field only. The corresponding transfer function (first two 
rows and columns of Eq. 1) is accordingly the conven- 
tional 2x2 Jones matrix [15]. In order to interpret the 
SOI via the conventional polarization parameters (diat- 
tenuation and retardance) , we then proceed to derive the 
Mueller matrix corresponding to this Jones matrix, us- 
ing standard relationship connecting Jones and Mueller- 
Jones matrices [15,16]. The resulting matrix is a diatten- 
uating retarder Mueller matrix characterized by diatten- 
uation (d(8)), retardance (5(9)) and orientation angle of 
the axes of the diattenuating retarder 
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M12 = dcos(2(/>),Mi3 = dsin(20) 
M 22 = cos 2 (2(f) + X COS (5) sin 2 (20) 
M23 = sin(20) cos(20) — xcos(5) sin(20) cos(20) 
M 24 = xsin(5) sin(20) 
M 33 = sin 2 (20) + xcos(S) cos 2 (20) 
M34 = x s'm(d) cos(20) , M44 = x cos(<5) , x = | ^1 — d 2 \ 

Here, the d and 5 parameters are related to E a (9), Ep(9), 
and thus to the amplitude scattering matrix elements, 
S 2 (9) and Si(9) as 
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In what follows, we analyze and interpret the competing 
three individual SOI effects via the pure Mueller matrices 
corresponding to each of these effects using the d and 
parameters encoded in its various elements. 



Case 1: d = 0,5 = ir <-> E a = or S 2 cos9 = 
—Si; Geometrical Berry phase effect 

The Mueller matrix (Eq. 3) corresponding to this special 
case represents pure geometrical Berry phase effect asso- 
ciated with azimuthal (0) rotation of polarization (spin) 
and subsequent generation of geometrical phase vortex 
(OAM) [1,7, 10]. The nature of SOI for such topologi- 
cal phase evolution can be understood by applying the 
Stokes vector corresponding to a horizontal linear polar- 
ization state (Si = [110 0] T ) on this matrix. The output 
Stokes vector becomes: S G = MS, = [1 cos 40 sin 40 0] T 
, which implies rotation of the incident linear polariza- 
tion state by twice the azimuth angle (20). This can 
be interpreted as, the two circular polarization modes 
(left and right), constituting the initial linear polariza- 
tion state, acquire opposite phase vortices (I = ±2) 
during their evolution in the scattering process [10, 14]. 
Any incident circular polarization state (e.g, LCP / 
RCP), on the other hand, undergoes flipping of helic- 
ity a = +1 <-> — 1 (evident from the negative sign of the 
matrix element M44 under this condition), and subse- 
quently generates phase vortex with topological charge 
I = ±2. This phenomenon has been observed in back- 
scattering from random medium [14]. Note, however, the 
condition S 2 cos9 = —Si can never be fulfilled for sin- 
gle scattering from dielectric Rayleigh scatterers (radius 
r « A, where scattering is primarily contributed by 
the lowest order TM scattering ai mode [17]) either for 
forward scattering (0 < 8 < tt/2) or for backscattering 
(ir/2 < 8 < tt) angles. In contrast, this condition can 
indeed be satisfied even for small angle (forward) scatte- 
ring from plasmon resonant metal nanorods [18]. In fact, 
all the competing SOI effects (discussed subsequently) 
can be significantly enhanced in controllable fashion, as 
illustrated below with selected examples. 

Briefly, metal nanorods exhibit two electric dipolar 
plasmon resonances, one at shorter wavelength (trans- 
verse resonance along the short axis) and the other at 
longer wavelength (longitudinal resonance along the long 
axis) [18]. In Figure 1, we show the wavelength variation 
( = 400nm — 800nm) of the computed d, and 5 param- 
eters for preferentially oriented silver nanorods having 
equal volume sphere radius r — 20nm and for two differ- 
ent aspect ratios of the rods (ratio of diameter to length) 
e = 0.65 and 0.95 respectively. The scattering angle is 
chosen to be 8 = 15° (as representative forward scat- 
tering small angle). The nanorods were oriented such 
that the long and the short axes of the rods are aligned 
along the laboratory X — Y (polarization) axes respec- 
tively [19]. The scattering matrices of the rods were cal- 
culated using T-matrix [19] method and the d and 5 pa- 
rameters were subsequently determined employing Eq. 
4. 

Several observations are at place. The diattenuation 
d parameter peaks at wavelengths corresponding to the 
two electric dipolar plasmon resonance bands (transverse 
at A ~ 425nm, longitudinal at A > 500nm) [18]. In con- 
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Fig. 1. (color online) Variations of diattenuation d (left 
axis, black lines) and retardance 5 (right axis, red lines) 
with wavelength A for silver nanorods having equal vol- 
ume sphere radius r — 20nm and for two aspect ratio 
values e = 0.65 (corresponding d and 5: solid lines) and 
0.95 (d and 6: dashed lines). The three different regions 
satisfying conditions for Case - 1, 2 and 3 SOI effects 
are marked as Region 1 (red), 2 (blue) and 3 (black) 
respectively. Conditions 1 and 2 are satisfied at two dif- 
ferent wavelengths for metal nanorods with e = 0.65 



(d^^t' and an d condition 3 is satisfied for metal 
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nanorods with e = 0.95 (d^ct an d ^mct)- The magnitudes 
of d (black dotted line) and 5 (red dotted line) for similar 
dielectric nanorods (r = 20nm,e = 0.65) are consider- 
ably weaker and do not exhibit spectral dependence. 
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trast, at the spectral overlap region of the two resonances 
(A ~ 450 — 500nm), d exhibits sharp concavity (d ~ 0). 
This rapid variation of diattenuation is a manifestation 
of preferential excitation of the two resonances with re- 
spective orthogonal linear polarizations for such prefer- 
entially oriented rods [20]. In addition to diattenuation, 
the nanorods also exhibit strong linear retardance S ef- 
fect, the magnitude of which peaks around the spectral 
overlap region of the two resonances. The enhanced 5 
effect can be attributed to the inherent phase retarda- 
tion between the two competing electric dipolar plasmon 
modes [20]. Importantly, the magnitudes of d and S can 
be desirably controlled by changing the aspect ratio e of 
the rods. In fact, at the overlap region of the two res- 
onances (A ~ 480nm), the parameter S and d can be 
tuned (by reducing e to ~ 0.65) so that they satisfy the 
condition 5 ~ n , d ~ (implying pure helicity reversal 
in forward scattering) (marked in red color, as region 1 in 
Figure 1). Clearly, such effect can never be achieved from 
dielectric Rayleigh scatterers (shown in Figure I), and is 
exclusively enabled by the presence of the two neigh- 
boring dipolar plasmon resonance modes in the metal 
nanorods. 



Case 2: d = 1, 6 ~ 0, <-> E a = ±E P or S 2 = 0/S 1 = 0; 
SAM to OAM conversion via pure diattenuation 
effect of scattering 

The Mueller matrix (Eq. 3) corresponding to this special 
case represents pure diattenuator matrix and signifies 
complete conversion of SAM to OAM (cr = ±1 -> 0, 1 = 
±1). This can be verified by applying the Stokes vector 
corresponding to input RCP state (Sj = [100 l] T ,cr = 
+1) on this matrix. The output Stokes vector repre- 
sents pure linear polarization state (cr = 0) : S D = 
[1 cos(2(f) sin(2<f>) 0] T . The generation of I = +1 vor- 
tex is manifested in the second and the third (linear 
polarization descriptor) elements of the output Stokes 
vector, by the appearance of the cos 2<j) and sin 2<f> fac- 
tors. The scattered light thus becomes completely lin- 
early polarized, carrying no SAM and accordingly the 
angular momentum is entirely carried by the OAM term 
(cr = +1 — y (7 = 0,/ = +1). As shown in Figure 1, 
the condition for such SOI (d = 1,6 ~ 0) can be satis- 
fied for plasmonic nanorods at wavelengths correspond- 
ing to the peaks of the dipolar plasmon resonances (at 
A ~ 525nm corresponding to the longitudinal dipolar 
resonance shown here, in blue color, as region 2). Note 
that such effect can arise even for scattering from dielec- 
tric Rayleigh scatterers [5]. The condition is fulfilled for 
Rayleigh scatterer at scattering angle 6 = 90° , where the 
amplitude scattering matrix element S2 vanishes, yield- 
ing d = I, 6 = 0. For large sized dielectric Mie scatterers 
(r > A) such complete conversion of SAM to OAM may 
occur at several other narrow range of angles depending 
upon the size parameter of scatterer [4,5]. In contrast, 
for the plasmonic nanorods, the phenomenon can be ob- 
served over a broad range of forward scattering angles, 
as illustrated in Figure 2. 

The angular variation of the derived d and S parame- 
ters for the silver nanorod with r — 20nm, e = 0.65, are 
shown here in Figure 2, for two different wavelengths, 
A 525nm (longitudinal resonance peak) and A 480nm 
(overlap region of two dipolar modes). Evidently, the 
condition d — 1,6 — 0(forA 525nm) is satisfied over al- 
most the entire range of forward scattering angles. 

Case 3: d = 0, 6 = 7r/2,o E a = iEp or 5 2 cos# = 
1S1; SAM to OAM conversion via pure retardance 
effect of scattering 

In this case, the resulting Mueller matrix assumes the 
form of a pure retarder matrix. For input horizontal 
linear polarization state (S, = [110 0] T ), the out- 
put Stokes vector of the scattered light becomes S G = 
[1 1 |(1 + cos 4</>) I sin Acj) sin 20] T , implying genera- 
tion of azimuthal angle <p - separated lobes of opposite 
circular polarization states. This effect is similar in na- 
ture to the Spin Hall effect of light, wherein an incident 
linear polarization state, evolve in different trajectories 
to generate spatially separated lobes of opposite circular 
polarization states (cr = ±1) [6]. Note in this case, the 
SAM to OAM conversion for input circular polarization 
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Fig. 2. (color online) Variations of diattenuation d (left 
axis, black circle) and retardance 5 (right axis, red cir- 
cle) as a function of scattering angle 9 for silver nanorods 
(r = 20nm, e = 0.65). The magnitude of d at A = 526nm 
(peak of longitudinal plasmon resonance) is ~ 1 for a 
broad range of 9 (satisfies condition for Case 2 SOI). In 
contrast, d for similar dielectric nanorods (r = 20nm, e = 
0.65) (black triangle) is quite low at smaller 9 and ap- 
proaches unity only at 6 = 90° . The magnitude of 6 for 
silver nanorod at A = 479nm (overlap region of the two 
dipolar plasmon resonances) remain close to 7r over a 
broad range of 9 (satisfies condition for Case 1 SOI). 

state (e.g., RCP, Sj = [1 1] T ), is similar to that ob- 
served for the previous effect. The output Stokes vector 
S = [1 - sin(2</>) cos(20) 0] T implies complete SAM 
to OAM conversion and subsequent generation of phase 
vortex / = +l(er = +1 — > a = 0, 1 = +1). From Figure 1, 
it is apparent that the condition for this type of SOI can 
be fulfilled in the spectral overlap region (A ~ 480nm) of 
the two dipolar plasmon resonances, for relatively larger 
aspect ratio of the silver nanorod (e = 0.95) (shown in 
blue color, as region 2 in figure 1). 

The results presented above show that enhancement 
of the SOI effects can be achieved by simultaneously ex- 
citing the two electric dipolar plasmon modes of metal 
nanorods. In general, such resonant enhancement of SOI 
should be possible in metal nanospheres also by exiting 
neighboring dipolar and quadrupolar plasmon modes. 
As an illustrative example, in Figure 3, we show en- 
hancement of the retardance S parameter by simulta- 
neous excitation of the electric dipolar (ai) and elec- 
tric quadrupolar (02) plasmon modes for a silver sphere 
with radius r — 50nm [17,18]. For this metal sphere, 
the ai and the ai plasmon modes are characterized 
by resonance peaks at 550nm and 365nm respectively. 
The results are thus shown for three different wave- 
lengths (A = 550nm, 365nm, peaks of the dipolar and the 
quadrupolar resonances respectively, and 400nm, over- 
lap region of the two modes). Apparently, 6 shows sig- 
nificant enhancement at A = 400nm, where both the 
a\ and a,2 plasmon modes are excited simultaneously. In 
contrast, a dielectric sphere having identical size (shown 
in the figure 3) does not exhibit any appreciable value of 
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Fig. 3. (color online) Variations of retardance as a func- 
tion of scattering angle 9 for silver nanosphere (r = 
50nm), for three different wavelengths, A = 365nm 
(quadrupolar resonance peak, black dashed line), 550nm 
(dipolar resonance peak, blue dotted line) and 400nm 
(overlap spectral region of the two modes, red solid line). 
Significant enhancement of magnitude of 5 is apparent 
at A = 400nm. The S- value at A = 400nm for a dielectric 
sphere having identical size (shown by circle) is rather 
low. 

this parameter, implying the corresponding SOI effect is 
rather weak in absence of the plasmon resonances. 

To summarize, the results demonstrate that each of 
the scattering mediated SOI effects can be resonantly 
enhanced by exploiting the interference of two neigh- 
boring plasmon modes in metal nanostructures (orthog- 
onal electric dipolar modes in rods or electric dipolar 
and quadrupolar modes in spheres). Importantly, the 
contributing SOI effects can be desirably tuned (opti- 
mized/enhanced) by changing the wavelength of light 
and controlling the size, shape of the nanoparticles. The 
developed generalized framework based on Mueller ma- 
trix approach enabled quantification and interpretation 
of each of the individual SOI effects, via their character- 
istic signature in the polarization patterns of the Mueller 
matrix elements (two-fold or four-fold azimuthal pat- 
terns depending upon the nature of the SOI), and via the 
diattenuation d and retardance S parameters encoded in 
the matrix elements. Although the present formalism is 
derived for scattering of plane waves, extension of this to 
include scattering of fundamental and higher order Gaus- 
sian beams, is also warranted. Since, the d and S param- 
eters can be directly determined from any experimental 
Mueller matrix [16], these parameters hold promise as 
novel experimental metrics for studying spin orbit inter- 
actions mediated by interactions of light with micro or 
nano-scale objects. 
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